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We present the possibility that Dirac and Majorana neutrino wave packets can be distinguished 
when subject to spin-gravity interaction while propagating through vacuum described by the Lense- 
Thirring metric. By adopting the techniques of gravitational phase and time-independent perturba- 
tion theory following the Brillouin-Wigner method, we generate spin-gravity matrix elements from 
a perturbation Hamiltonian and show that this distinction is easily reflected in well-defined gravita- 
tional corrections to the neutrino oscillation length for a two-flavour system. Explicit examples are 
\£) • presented using the Sun and SN1987A as the gravitational sources for the Lense-Thirring metric. 

f"^ | This approach offers the possibility to determine the absolute neutrino masses by this method and 

identify a theoretical upper bound for the absolute neutrino mass difference, where the distinctions 
between the Dirac and Majorana cases are evident. We discuss the relevance of this analysis to 
the upcoming attempts to measure the properties of low-energy neutrinos by SNO and other solar 
neutrino observatories. 
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> 1 I. INTRODUCTION 

co . 

As arguably the most elusive of the known subatomic particles in existence, neutrinos nontheless offer great promise 
\Q i for revealing insights at both the microscopic and cosmological level. Because they are generated solely from the weak 
interaction and have an extremely small scattering cross section, to just even observe a neutrino event is a major 
experimental accomplishment. Neutrinos have the surprising property of having exclusively negative helicity due to 
parity violation of the weak interaction, and that only their left-handed chiral projections interact with matter. For 
all these reasons, it was understandable to assume that neutrinos were strictly massless particles like photons which 
propagate at the speed of light, and are represented as such in the Standard Model of particle physics. However, there is 
| | clear evidence that the Standard Model is incomplete because of the recent experimental discovery by Superkamiokande 
[1] and the Sudbury Neutrino Observatory (SNO) [2] demonstrating that neutrinos can transform from one species to 
another while propagating through space. This property, known as neutrino flavour oscillation [3], strongly suggests 
that neutrinos have non-zero rest mass, since the standard treatment of the theory for a two flavour system requires 
that each flavour has a distinct mass eigenstate. 

The current theory of neutrino oscillations is based on solving the Dirac equation in flat space-time for the mass 
eigenstates and their associated energy eigenvalues. One immediate consequence is that we lack sufficient information 
to determine their absolute masses outright, since only the mass-squared difference can be determined directly by 
observation. Because of this fact, it is not obvious how to obtain the absolute neutrino masses without resorting 
to debatable theoretical speculations. Since an enormous number of neutrinos were produced in the early Universe, 
having a reliable figure for the neutrino masses can give cosmology researchers better insight about the nature of dark 
matter, with implications for the Universe's overall evolution. There is also a conceptual difficulty with the current 
approach to neutrino oscillations, in that neutrinos are largely produced in supernovae and within the stellar cores 
of astrophysical sources, where for the case of a neutron star the curvature of space-time near its surface is close to 
the Schwarzschild limit. Therefore, a more complete approach should be within the framework of a curved space-time 
background, as described by Einstein's general theory of relativity. 

On a different level entirely, neutrinos can be modelled as the propagation of plane waves, wave packets, or quantum 
fields. The obvious advantage for using plane waves comes from the simplicity in performing calculations, but at the 
expense of losing possibly significant physical insights. Indeed, the fact that neutrinos are quantum objects which 
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can be localized in a detector suggests that they should at least be treated as wave packets [4] with some spatial 
extension subject to the Heisenberg uncertainty relations. This point becomes especially relevant when considering 
the effects of gravitation on neutrino wave packets, since both are inherently non-local in nature, where space-time 
curvature may have a non-trivial influence on their intrinsic properties as they propagate through space. Conversely, 
it is theoretically possible to glean insights about the nature of space-time at the quantum mechanical level that may 
provide some observational clues towards attaining a viable quantum gravity theory. In fact, much research has been 
done to better understand the interface between quantum mechanics and classical gravitation, where most recently 
it was shown [5] that the gravitational field can distinguish between the helicity and chirality of massless spin- 1/2 
particles. This is in contrast to the prevalent understanding [6] that the helicity and chirality operators coincide for 
massless spin- 1/2 particles in a flat space-time background. 

Another more fundamental debate concerns whether neutrinos exist as Dirac particles, where the antineutrinos 
are physically unique particles, or Majorana particles, where the neutrinos and antineutrinos are identical. More 
precisely, this means that Dirac neutrinos and antineutrinos have a distinct set of quantum numbers, while the 
Majorana set share the same quantum numbers. While lepton number is a conserved quantity for Dirac neutrinos, it 
is violated for Majorana neutrinos. Currently the only proposed method to distinguish between the two types comes 
from neutrinoless double beta decay [7, 8], which can only happen for Majorana neutrinos. Because of the extreme 
difficulty in making such an observation, it may be a long time before this ambiguity can become clarified by this 
approach. As for neutrino oscillations, current understanding suggests that it is not possible to distinguish between 
Dirac and Majorana neutrinos by this method. This is because both neutrino types share the same left-handed 
projection of their respective four-spinors [9], which generate the same physical predictions in this context. Again, 
however, it is important to emphasize that this understanding comes about in the absence of a gravitational field. 

In this paper, we consider the possibility of whether classical gravitation, described in terms of the Lense-Thirring 
metric [10] can have an impact on neutrino wave packet propagation, with consequence for neutrino oscillations. 
This is done by treating the gravitational field as a phase contribution to the neutrino wavefunction leading to an 
interaction Hamiltonian. By using the techniques of time-independent perturbation theory, we can determine the 
perturbation energy from the matrix elements corresponding to the propagation of neutrino wave packets, which 
are represented by a power series expansion in terms of their rest masses. A preliminary investigation [12] into this 
problem following a semi-analytical approach shows that it is theoretically possible to extract the absolute neutrino 
masses from direct observation because the oscillation length becomes dependent on both the linear mass difference 
coupled to the gravitational field, as well as the mass-squared difference. 

Since this first investigation considered the problem exclusively for Dirac neutrinos, it left open the question of 
whether Majorana neutrinos would behave differently from their Dirac counterparts under the same conditions. The 
purpose of this paper is to address this question by performing the same general analysis as before, but with much 
greater attention to detail emphasized throughout. As a result, we have a set of calculations to present with virtually 
no approximations employed, in which we show that the gravitational field can fundamentally distinguish between 
Dirac and Majorana neutrinos [13]. We begin in Sec. II with the formalism on neutrinos in curved space-time, the 
Lense-Thirring metric, and a description of gravitational phase required for the time-independent perturbation theory 
calculations in this paper. This is followed by Sec. Ill, which provides both the formalism of neutrino wave packets. 
The resulting matrix element calculations are then described by Sec. IV for both Dirac and Majorana neutrinos. After 
this, we proceed in Sec. V to determine the mass-dependent perturbation energy via the Brillouin-Wigncr approach 
normally used in condensed matter theory. This has the advantage of generating exact energy corrections for a two- 
level spin system, eventually leading to the main formal expressions in this paper. We then perform some extensive 
numerical analysis in Sec. VI using both the Sun and the SN1987A supernova as examples to predict the gravitational 
corrections on the neutrino oscillation length, and present a method to extract the absolute neutrino masses by this 
approach. Finally, we conclude in Sec. VII with a discussion about possible future developments that may follow from 
this investigation. 

II. NEUTRINOS IN CURVED SPACE-TIME WITH GRAVITATIONAL PHASE 

A. Covariant Dirac Equation 

For all calculations to follow, the space-time metric has —2 signature and we adopt geometric units G = c = 1 [14], 
such that all physical quantities are described in units of length defined by M, the mass of the gravitational source. 
We start with the covariant Dirac equation in curved space-time [15] for spin-1/2 particles with mass m: 

il>{x) = 0, (1) 



i^{x)D ll - 
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where the set of space-time dependent gamma matrices {7 M (x)} satisfy {7^(2;), 7^ (x)} — 2 g M!y (x), and the covariant 
derivative operator = + iV ^ is described in terms of spin connection F^, such that J v (x) = 0. By defining 
a local Minkowski frame at each space-time point on the manifold for a given set of orthonormal tetrads {e^} and 
basis one-forms {e' 1 } labelled by hatted indices and satisfying (e^,ej>) = S^o , the curved space-time metric is 
g = f)iiv ef 1 <8> e u , where r)p,p is the Minkowski metric tensor. The general metric tensor g^ v is then related to r)p,p in 
terms of vierbein projections {e"^} , {e^a} satisfying e a — e a p e 13 and e& = e^a ep, such that 

6 ^ $ — ^ /3' ^ ot 6 v — ' u , (2) 

9^ = r n & ^e & ll eK- (3) 

The spin connection matrix is then 

r„ = -\^{x)Y a ^ = -I^r^e^, (4) 

where a a/3 = f [7", 7^] are the Minkowski space-time spin matrices, and are the Ricci rotation coefficients derived 
from the Cartan equation of differential forms 

de^ + r^.e^Ae' 5 = 0. (5) 



B. Dirac Hamiltonian in Lense-Thirring Space-Time 

The Lense-Thirring metric in cartesian co-ordinates x M = x 1 , x 2 , x 3 ) = (t, x, y, z) is given by 



1 - 



2M 



r J 
4 MttR 2 



dt®dt- 



1 + 



2M 



r) 



(dx <S) dx + dy <£> dy + dz <g> dz) 



+ 7 



[a; (dy <g> dt + dt <x> dy) - y (dx (g) di + dt <g> da;)] 



(6) 



used to describe a weak gravitational field in the vicinity of a rotating star, where R is the radius of the gravitational 
source, f2 is its rotational velocity, and r = \J x 2 + y 2 + z 2 . Given that a and (3 are the usual four-dimensional 
Dirac matrices and 53 J = (T 0j is the spin angular momentum operator, the Lense-Thirring spin connection components 
following (4) are then 
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Therefore, the corresponding Dirac Hamiltonian derived from (1) up to leading order in M/r is 
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(11) 



where L 2 = xp$ — yp x is the orbital angular momentum operator in the z-direction. For future reference, (11) is 
regarded as the zeroth-order part of the total Hamiltonian with an interaction term, where the latter is a perturbation 
due to the gravitational phase contribution, whose details are listed immediately below. 
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C. Gravitational Phase and the Interaction Hamiltonian 



For this investigation, we want to consider the effect of gravitation on the quantum propagation of a neutrino 
wave packet from the gravitational source to the Earth. Such a condition implies that, unlike other scmiclassical 
treatments often found in the literature [16], we need to average over all possible trajectories within phase space that 
are accessible by the neutrino, as opposed to integrating over its classical trajectory described by the geodesic equation. 
Such a consideration requires that we use time-independent perturbation theory, where the interaction Hamiltonian 
is dependent on the space-time metric. This can come about following the introduction of the gravitational field via a 
phase shift in the wavefunction. The gravitational phase <£>g for a weak gravitational field is then described by [17-20] 



<&G = \ I dz x h Xa (z)p a - - [ dz x [h Xa . fj (z) - h xf} . a (z)] L af3 
2 JxZ 4 Jxt! 



(12) 



where = g^ v — rj^ is the metric deviation, p^ and L a>3 are the momentum and orbital angular momentum 
operators of the free particle, and z^ = (t', x' , y', z') is the integration variable over some arbitrary classical path in 
space-time. One of the main properties of the gravitational phase in this form is that a wavefunction described with 
(12) becomes a solution of the weak-field covariant Klein-Gordon equation, exact to first order in hp V and invariant 
under h^ u — > h^ v — (^, v + the co-ordinate gauge conditions which follow from the co-ordinate transformations 

When applied to a closed space-time integration path, (12) leads to the covariant generalization of Berry's phase 
[21] with useful applications for calculations involving particle intcrfcrometry [22]. This form of the gravitational 
phase is comparable with other formulations [23-26] available in the literature. It has also been shown [27] that use 
of (12) correctly reproduces the degree of particle deflection predicted by general relativity in the geometrical optics 
limit. When $g is applied to spin-1/2 particles accelerated in storage rings for a metric describing noninertial motion, 
it gives rise to all gravitational-inertial effects that have been either observed outright [28, 29] or are very likely to 
exist [30, 31]. 

For the purpose of this investigation, the gravitational phase can be introduced by means of the transformation 
ip(x) — > exp (i$c/ft) ip{ x ), leading to the total Hamiltonian H = H + -ff$ G , where 



H$ G = a • (V$ G ) + (V t $ G ) 



(13) 



is a first-order correction to Hq. Therefore, (13) is the interaction Hamiltonian sought after, which describes a direct 
spin coupling to the gravitational field via $g , where the gravitational phase can be written as 
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(14) 



In terms of the Lense-Thirring metric, the explicit expressions for V^c&g are 
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The orientation of the neutrino beam in space, expressed in terms of spherical polar co-ordinates, is shown in Figure 1. 

While this treatment does not pose any problems for the spatial separation r — r' found in (16)— (18) , the time 
interval t — t' requires more careful consideration, given that the spatial part of the event x^ = (t, r) is not constrained 
to satisfy a classical trajectory defined by the Einstein field equations. However, there is a way to at least approximate 
t — t' quantum mechanically by imposing Poincare invariance on Ax^ = x^ — x'^ for a timelike interval At, such that 



(At) 2 = (At) 2 - (Ar • Ar) w (At) 2 



(20) 



where 7 = E/m, and E — \/(p ■ p) + m 2 . It follows that 
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(21) 



Even with this approximation, however, there is a technical challenge in performing the integration over all spatial 
angles contained in (21), so we further approximate this time interval by assuming 
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where 



(Ar ■ Ar) = j (r - r') ■ (r - r') dfl = 4tt (r 2 + r' 2 ) 



is the angular average over dfi = sin^d^dtp. This construction leads to the expression 
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(24) 



which we substitute into (15)-(18) for further evaluation. It immediately follows that (24) identically vanishes in the 
limit as x'^ — > x^. 



III. NEUTRINO WAVE PACKETS AND SPIN-GRAVITY INTERACTIONS 

There is some debate [4] over whether neutrinos exist strictly as plane waves with precisely defined momentum 
eigenstates, or as a linear superposition of plane waves with some momentum spread. While the computational 
simplicity of the plane wave approach over using wave packets is obvious, the problem with such an assumption is 
that the neutrino is not localizable in space, according to the Heisenberg uncertainty principle. Although the choice 
to use wave packets is motivated more by the need to see the effects of curved space-time on some localized region 
where the neutrino can be identified, we acknowledge that it makes intuitive sense to treat neutrinos like any other 
quantum particle according to standard quantum mechanics. Therefore, we assume that the neutrino wavefunction 
is a Gaussian wave packet composed of a linear superposition of plane waves in momentum space, described by 



(2tt) 3 / 2 



Jd 3 kak)e lk ' r \U(k)), (25) 



where \U(k)) is a free-particle spinor with momentum p = hk, and the normalized Gaussian function of width a p 
and mean momentum (p) — hk is 



m = V^fT 2 6XP 



(fc-fco) 2 
4a 2 



(26) 



In order to make use of time-independent perturbation theory consistently, we need to assume that tx p is constant. 
Of course, we expect that a realistic neutrino wave packet will spread with time, but assume in this model that the 
spreading is so rapid that it will not adversely affect the main results to follow. Therefore, we assume implicitly 
that the perturbation energy emerges over a very short time scale and remains frozen almost immediately after the 
neutrinos start propagating through space, when the gravitational potential rapidly decays to zero. 
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FIG. 1: Orientation of the neutrino beam in relation to a rotating gravitational source of mass M, radius R, and angular 
velocity f2 which describes the Lense-Thirring metric. As a quantum mechanical particle, the neutrino has the dynamical 
freedom to occupy any position r' — (r , 9' , tp') while propagating through space from the source to the detector defined by the 
co-ordinates r = (r,9,tp). 



The matrix element for a transition involving neutrino wave packets due to spin-gravity coupling is then formally 
described by 

(^(r)\H^(r)) = ^3 j dV d 3 fe d 3 fe' £(k) Z(k') exp[i (k - k') ■ r'] (U(k')\H* G (r,r')\U(k)), (27) 

where we integrate over all phase space, excluding the region occupied by the gravitational source. It is from the 
choice of \U(k)) in (27) which determines the nature of the matrix element due to spin-gravity coupling, and hence 
the degree of mass-dependent energy difference predicted for the Dirac and Majorana neutrino. This comparison 
between Dirac and Majorana wave packet constructions and their impact on (27) immediately follows below. 



A. Dirac Neutrinos 



Assuming the Dirac representation [6] for a: and /3, the corresponding free-particle Dirac spinor for use in (25) is 



\U(k))f^ = sj^^ ^ j ® |±), (28) 

where 

cos (6/2) \ / -sin (9/2) \ 

sm(e/2)e l(p J ' 1 ; \cos(9/2)e^ J { > 

are the two-component spinors [11] with positive and negative helicity, respectively, whose spin quantization axis is 
defined in terms of (19) along the neutrino beam's radial direction. Then the interaction Hamiltonian (13) in the 
Dirac representation is 

H^ G = (a ■ V$ G ) ® ( \ \ \+ (V t * G ) ® ( I J ) ■ (30) 
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It is straightforward to show that the transition from \U(p)) to \U(p')) to evaluate (27) is described by 
(U(p')\H iG \U(p)) m - = 
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(32) 



is the gradient of <&g symmetrized with respect to p and p' for the purpose of integration over all intermediate 
momentum states. While this symmetrization of V m $q nas n0 bearing on the final evaluation of (27), it does serve to 
help identify in (31) the symmetric and antisymmetric contributions of the final matrix clement under the interchange 
of p and p'. We note immediately that the symmetric part of (31) is independent of the initial spin polarization of the 
neutrino, while the antisymmetric part is coupled to the spin operator. This highlights an important feature about 
the wave packet treatment, in that a spread of momentum states in the neutrino allows for the possibility to register 
a non-trivial contribution of the Dirac matrix clement due to spin flip. In contrast, it is impossible for a plane-wave 
neutrino to see the same effect, since the spin flip part of (31) identically vanishes as p — > p'. When integrating over 
phase space, it follows from (31) that the Dirac matrix element is 



Mr)|ff* G Mr)> 
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, (33) 



where the spin flip transition corresponding to the neutrino beam's orientation in space is represented in terms of a 
cartesian co-ordinate frame by 



(±|<t|=F) = [cos 9 cos (p ± i sin ip] x + [cos 9 sin ip + i cost/?] y — sin 6 z. 



(34) 



B. Majorana Neutrinos 



Performing the same set of calculations for the Majorana neutrino requires some preparation beforehand. This is 
because the Majorana equivalent of the frec-particle spinor (28) must be self-conjugate under the charge conjugation 
operation 



\u c (k)) ± = c 7 6 |tr(fc)}±, c = z 7 V- 

For a Majorana spinor \x), this implies that [8] 

lx c ) = ±lx>- 



(35) 
(36) 



In order to best satisfy this requirement, we need to reformulate the problem in terms of the Wcyl or chiral represen- 
tation [6] for the gamma matrices. The reason follows from knowing that the Dirac spinor (28) can be decomposed 
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into the sum of a left- and right-handed projection expressed as more primitive two-component spinors. Then the 
Majorana spinor can be constructed as the sum of the Dirac spinor's left-handed projection and its charge conjugate 
under (35), which then satisfies (36) automatically. 

We can explicitly proceed by first making use of the unitary matrix U which relates the Dirac and Weyl represen- 
tations according to 

7{kyi = U^ iiac U\ (37) 

where 



"-M'-^-M-iD (38 » 



in the Dirac representation. With (38), the interaction Hamiltonian in the Weyl representation becomes 

H* G = (<x ■ V$ G ) ® ( J ° ) + (Vt$ G ) ® ( I ? V (39) 



o -l j ' y ^ y o l / 

A simultaneous transformation of (28) then leads to Weyl form of the free-particle Dirac spinor 



where the left- and right-handed two-component spinors are 
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W ± - |(1 + ^ 5 ) HWff" - £ ( (' + j^) ) * l±>, (41b) 

such that |[/(fc))^ cyl = (\vl) + \vr))±- By then applying (35) to \vi) and \vr), we obtain the charge conjugated 
spinors 

- ^^(-,V+^)) Sl±> = 71^(-(i-k)^) 8l±> ' <42b) 

where (cr?) = — er* er^ er^. Therefore, we can define the free-particle Majorana spinors for the matrix element (27) as 

114^)}^ =(K> + K»±, (43a) 
\W 2 (k))™^ ={\u R )-\u%)) ± , (43b) 

which by construction automatically satisfy (36). That is, 



V2V 2E y -(l-^) J " y/2 V ^ _ (l _ » 



(44a) 



hcr k\ \ . /-= / i ft<xfc 



The most transparent way to evaluate the Majorana matrix clement for (44a) and (44b) is to determine each of 
the components using (41a)-(42b). For transitions from a left-handed to left-handed (right-handed to right-handed) 



9 



two-component spinor, we obtain 
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where the upper (lower) sign corresponds to the left- (right)-handed two-component spinor. We immediately verify 
from comparison with (31) that 



(v l ( P ')\Hz g \v l (p)) + (Mp')\H* g \Mp)) - (U(p')\H* G \U(p)} 



Dirac 



(46) 



as expected. Similarly, the conjugate left-handed to conjugate left-handed (conjugate right-handed to conjugate 
right-handed) two-component spinor transition is 
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+ TO 



£7 + TO 



(V t * G ) s + 



£ + TO 



+ i cr 



(V t $ G ) s (p x p') /((V$ G ) s xp) ((V$ G ) s xp') 



(£ + to)(£' + to) 



£ + TO 



E' + m 



±i 



. (V$ G ) S • (p x p') 



(-E + m)(£' + m) 



(47) 



Therefore, from using (45) and (47), it follows that the relevant transition matrix element for the Majorana neutrino 
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is 



(W m (p')\H* G \W m (p)) 



Maj. 



(U(p')\H* G \U(p)) 



Dirac 



2E 



I E + m E' + m 
+ \l \l — Off' 



2£' 



(V t $ G ) g (p x p') / ((V* G )s x P) ((V*g) s x P') 



(E + m)(S' +m) 



£ + m 



E' + m 



(V t $ G ) s (p x _ / ((V$ G ) S x pf _ ((V<P G ) s xp'f 
(E + m) (E' + m) y E + m E' + m 



(p-p') 



(E + m)(E' + m) 



[(<T-V$ G ) S -^(V,$ G ) S ] 



+ ( (V^c) s + ^k^l) i'-P-***) + ( {Vt , G)s + ((^ok -P)\ {cr-p'-a* p 



E' + m J E + m \ E + m 

and when integrated over all phase space, the Majorana matrix element is then 

«>i( 2 )(r)|!^ G |^i( 2 ) (r)} Ma j- = Mr)|ff* G Mr)> 



"A 



E' + m 



1 /'.iVuW^' £(fc)£(fc') /i? + m /s ' + m 



(2^)3 



2E V 2£ 



(48) 



x \ sin [(fc - fc') ■ r'] 



(±|ff|T) 



ft 2 (V t $ G ) s (fc * fcQ _ / fe((V$G) s xfc) ft((V$ G ) s xfc') 
(-E + to)(£' + to) I £ + to + m 



(±|ff|T) S 



^ (V f $ G ) s (fc x fc') fi 
(E + m)(E' + m) 



'fi((V$ G ) g x fcf fi((V^ G ) s x fc 



n9 



£ + m 



£' + m 



± cos [(fc - fc') • r'] (±| 



ft 2 (fc • fc') 



(£■ + m)(E' + to) 



[(<r.V$ G ) s -ff 5 (V„<!>G) E 



+ hi (V 4 $ G )s + 



. fe ') \ (er • fc - cry k y 



E' + m 



E + m 



M(V t #G)s 



h((V^ G ) s -k) 
E + m 



cr ■ fc 



rVk' y 



E' 



(49) 



Immediately, we can recognize that the Majorana matrix element (49) departs significantly in form compared to 
its Dirac counterpart (33). Although a more detailed analysis follows later in this paper, it is already clear that the 
outcome of using a self-conjugate spinor is to induce a frame-dependent expression which favours the y-direction. 
This comes from the presence of the cr v Pauli matrices in |f£) and \iy%), which contribute exclusively to (47). Indeed, 
if we tag cr y in (42a) and (42b) with a parameter e, then (47) is of order e 2 and vanishes in the limit e — > 0, leaving 
only (45) to contribute to (49), as expected. 

The Dirac matrix element (33) and the corresponding Majorana version (49) for iVt^i)^ 1 ^' are the two physically 
relevant cases that we consider exclusively for the rest of this paper. However, for the sake of completeness, we also 
briefly examined the matrix element for a transition from |Wi)± aj ' to |W2)± aj '. To do this, we need the transition 
from the conjugate right-handed two-component spinor to the left-handed (conjugate left-handed to right-handed) 
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expression, where we obtain 

^L(R){p')\H^ G \v c R(L) {p)) 



1 E + m E' + m 



2 V 2E V 2E 



(<r x V$ G )i-i(V y $ G ) s 



1 - 



(p-p') 



y (V^G) S ■ (P X P ') 



(E + m)(E' + m) J (E + m)(E' + m) 



+ 



(<r x p) y - ip$ 



E + m, 



( v,* o)s+ «^,- 



(er x p'f -ip ,?} 



(V t $ G ) s + 



((V^g) s -p) 

£ + m 



± 



(V t $ G ) s (pxp'f /((V* G ) s xp) fl ((V$ G ) s xp'f 



(£ + m)(i;' + ra) I E + m 



E' + m 



+ i 



(V 4 $ G ) S [o- x (p x p')f [ffx ((V* G )s x P)]" [o- x ((V*g) s x p')P 



(E + m)(E' + m) 



E + m 



E' + m 



+ ia y 



(V t $ G ) s 1 



(p-p') 



f ((V$ G ) s -p) , gV^gisV) 



{E + m){E' + m)J \ E + m E' + m 

Therefore, it follows that the Majorana transition matrix element from (44a) to (44b) is 
(W 2 {p')\H^ G \W 1 (p)) M ^ = 



(50) 



— i 



. E + m E' + m 



2E V 2E 



(V t $ G ) s [<x x (p x p')f 
(E + m)(E> +m) 



[cr x ((V*g) s x P)] v [tr x ((V$g) s x p' 



£ + TO 



E + m E' + m 



2E V 2£ 



(V„$g) s 1 



(p-pO 



(£ + m)(£' + m) 



+ (V t * G )s + 



«V$gW)\ p 



£' + m J E + m 



(V*$g) s 



((V$ G ) s -p)\ p 



£ + m E' + m 



E + m E' + m 



2E V 2£ 



(VtjG)s (P x pT ((WcWf ((V^g) s x p') 



(E + m)(E' + m) 



E + m 



E' + m 



(51) 



What is most surprising about this result is that a real and non-trivial contribution survives due to its first three 
terms in (51). This is because they are both antisymmetric under the interchange of p and p' with an overall factor 
of i, which when multiplied by exp [i (fc — fc') ■ r'] and integrated over all phase space, yields a real contribution that 
comes from sin [(fc — fc') • r'] in the integrand, as demonstrated in (49). 
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IV. EVALUATION OF THE GRAVITATIONAL PHASE MATRIX ELEMENT 

A. Formalism 

To evaluate the Dirac and Majorana matrix elements (33) and (49), respectively, we make use of spherical symmetry 
to integrate over phase space in terms of spherical polar co-ordinates. Then the wave number k is expressed in the 
standard form 



k x = k sin A cos^i, k y = k sin A sin/z, k z = fccosA, 



(52) 



where A and fi cover the unit sphere and k goes from zero to infinity. Since we want to determine gravitational 
corrections to the neutrino oscillation length in terms of neutrino mass, we proceed to expand (33) and (49) by power 
series in m, up to second order. The next step is to integrate over all the position and momentum space angles. 
From a first impression, this seems like a formidable challenge considering that the plane wave term exp [i k ■ r'] is 
angle-dependent. Fortunately, this problem can be circumvented by employing the Raylcigh plane wave expansion 
[32] 



jk-r' 



47r E E i l 3i(kr')Y^(X,^ Y lm {6, V ) 



(53) 



1=0 m= 



where ji(kr') are the spherical Bessel functions and Y\ m are the spherical harmonics for angles in position and 
momentum space. Since the angles in (33) and (49) are expressed as polynomials of sinusoidal functions, they can be 
expressed as a series expansion of spherical harmonics to take advantage of their orthogonality relations. It follows 
that the series expansions truncate by virtue of the orthogonality conditions and the integrations over angles are 
performed exactly for both the Dirac and Majorana matrix elements. 

For the final integrations to evaluate, we first proceed to integrate over k, in which the integrals take the general 
form 



Qc(n) = / 
Jo 

Qs(n) = f 
Jo 



exp 



exp 



(fc-fcp) 2 



cos(fcr')fc™ dfc, 
sin(fcr')fc" dfc, 



(54) 
(55) 



where fco = |&o| an d n is an integer with — 2 < n < 2. All details related to evaluation of (54) and (55) are contained 
in Appendix A, while the non-trivial radial integrals which remain are listed in Appendix B. For future reference, the 
width of the wave packet is expressed in terms of the dimensionless parameter q = k n /a p , and the neutrino mass is 
now in units of its mean momentum, such that we have the dimensionless quantity m = m/(hk n ) <C 1. 



B. Dirac and Majorana Matrix Elements 

For evaluating (33), we obtain the expression 

(^(r)|ff $G |^(r)) Dirac = (ftfcb){y [Co + Cr m + C 2 m 2 ] + sinfl [fl + flim + D 2 m 2 ] J 



for the Dirac matrix clement, where 



/ C,(r')dr', Dj = / D,(r')dr' 

JR J R 



(56) 



(57) 
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are dimensionless functions dependent on q, ko, R, and r. The integrands for (57) in terms of (54) and (55) arc 
Co(r') = yj^ 3 |-^ 7 [Q c (l)Q c (2) + g s (l)Q s (2)] 

+ T~i2 1 3 0c(2) + 5 Qc(i) 0,(i)] - t^73 0,(o) 0,(i)| . ( 58a ) 

= \ {^{£) {? [Qc(l) 2 -Qs(l) 2 + 3Qc(0)Q c (2)-3Q s (0)Q s (2)] 

- [3 Q c (2) Q 8 (-l) + 5 Q c (0) g s (l)] + ^ [5 Q 8 (-l) 0,(1) - Q s (0) 2 ] } , (58b) 
C2(r') = ^/J(^) {^[3Q c (-l)O c (2)-7O c (0)Q c (l) + 3O s (-l)O s (2)-7Q s (0)O s (l)] 

- I 3 Qc(0) 0,(0) - 5 Q c (-1) 0,(1) - 3 Q c (2) Q s (-2) + 5 Q c (l) Q a (-1)] J . (58c) 

A,(rO = ^yj(±) 8 {^[O c (l) 0.(2) + Q C (2) 0.(1)] 

4 r 2 1 fi r 2 8 r 2 1 

- 7T^3 [ 2 ^c(l) 2 + 3 0c(0) Q 5 (2) + 3 Q s (l) 2 ] + —j Oc(l) 0.(0) - O,(0) 2 , (58d) 

IvO T fvQ T Kq T J 

Hr') = ^/5(j) 3 {-^[Oc(0)O s (2)-Oc(2)O,(0)] 

+ IS" t 4< W°) + 3 0-(-l) 0,(2) - 3 0,(0) 0,(1)] 

- 7? [0c(o) 0,(0) + Oc(i) 0-(-i)] + 04-1) 0,(0)} , (58e) 
^ 2(/) = \/J (£) {~^~ [<3c(_1) 0s(2) " Oc(0) 0s(1) " 0c(1) Os(0) + 0c(2) 0s( ~ 1)] 

+ -^g- [4 Oc(-l) Oc(l) - 4 Oc(0) 2 + 3 0»(-2) 0,(2) - 3 Q,(0) 2 ] 
4 /rn r 2 

- -^r- [Qc(-i) 0,(0) - 2 Oc(0) 0-(-i) + Oc(i) 0»(-2)] 

+ ^|^[O,(-2)O,(0)-O s (-i) 2 ]}. (58f) 

Substitution of (54) and (55) into (58a)-(58f) followed by the radial integration then leads to explicit expressions of 
Cj and Dj, which are found in Appendix C. 

Similarly for the Majorana matrix clement, we obtain from (49) the expression 

(4>i(2) {r)\H* a \4> m (r)) Ma J- - {hko){y [C + d m + C 2 m 2 ] 

± sine sin^ — (±|<t|t) S [C y + C ly m + C 2 g m 2 ] + — ^ — <±|«t|=f>* [Dost + D ljb m + D 2£ fh 2 ] \, 

(59) 

where 



0?5 



poo roo 

/ C y{r')dr', D j& = / D j£ (r')dr' (60) 
J r Jr. 
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and 



Coy(r') = ^ (£) {^[Qc(l)Q c (2) + Q s (l)Q s (2)] 

- K C ( 2 ) Qs(0) - Qc(l) 0,(1)] - ^3 0,(0) 0,(1)} , (61a) 

(5ls(r ' ) = 6^ (j^) {7 t^^ 1 ) 2 ~ Q^ 1 ) 2 ~ Qc(0) <9c(2) + 0,(0) 0,(2)] 

+ ^ [Q c (2) Q 8 (-l) - Q c (0) Q s (l)} + ^ [Qs(-l) 0,(1) - Qs(0) 2 ] } , (61b) 
C**^) = {-^[ZQMQc{i) + ZQ s (V)Q s (i) + Qc(-i)Qc{2) + Q s {-i)Q s {2)] 

+ t-Qc(-l) 0,(1) - Qc(0) 0,(0) + Q c (l) Q a (-1) + Q c (2) Q.(-2)] 

+ [<W" 2 ) + 3 Q'(-V Q-(°)]} . ( 61c ) 

A»(rO - |/5(^) 3 {^3^(i)Oc(2)-g 5 (i)o s (2)] 

- ^4 [Qc(2)Q s (0) + 4 0c(l) 0,(1)] + £^5 0,(0) 0,(1)}, (61d) 
^ (r,) = iV?(^) 3 {~^ [Oc(0)Oc(2) + 3Oc(l) 2 + 0,(0) 0,(2) ~0,(1) 2 ] 

+ I 4 Qc(o) Q'O-) + 3 Qc(i) Qs(o) + q c (2) 04-1)] 

- l^Q s (-l)Q s (l)J, ( 6 le) 

^2x(r') - (£) {^[-Oc(-l)Oc(2) + 13Oc(0)Oc(l) + O s (-l)O s (2) + 3O 5 (0)O,(l)] 

- -75- h 4 0c(-i) 0,(1) + 70 c (0) 0,(0) + I00c(i) O s (-i) - Oc(2) Q.(-2)] 

- i^l [O s (-2) 0,(1) - O s (-i) 0,(0)] J . (6if) 

The final expressions for Cj y and Dj£ are also listed in Appendix C. 

Even at the matrix element level, we notice that the expressions (56) and (59) for the Dirac and Majorana neutrinos, 
respectively, have fundamental properties worth noting. Regarding the Dirac matrix element, we identify the Cj and 
Dj in (56) as the spin-diagonal and spin-flip structure functions for j = 0, 1, 2, since Cj is coupled to M/r and Dj is 
coupled to MilR 2 /r 2 . The fact that sin^ is also coupled to the spin- flip terms clearly indicates that only terms with 
the ^-component of (34) contribute to the matrix clement. This leads to an obvious interpretation that rotational 
inertia from the Lense-Thirring metric induces the helicity to flip for any neutrino propagation that is not along 
the source's axis of symmetry. The effect is analogous to that of an inhomogeneous magnetic field which forces the 
particle's spin to flip. Perhaps most importantly of all, (56) has terms which are linear in m, which come from 
the normalization coefficient y/ (E + m)/(2E) in |[/(fc)) Dlrac . This property becomes most relevant when considering 
gravitational corrections to neutrino oscillations. 

As for the Majorana matrix element (59), it shares the same type of properties as found in (56), but with obvious 
differences. First, while the spin-diagonal terms Cj are common to both matrix elements, the remaining terms have an 
overall factor of sin 9 sin ip, corresponding to the y-component of the neutrino beam. This is a direct consequence of the 
Majorana neutrino's self-conjugation condition, since the charge conjugation operation has er^ present in its definition 
[8], and suggests a preferred direction orthogonal to the source's symmetry axis. However, this is inconsistent with 
the fact that the Lense-Thirring metric is axisymmetric, and implies that the 93-dcpcndence on (59) is purely due to 
how we defined the co-ordinate axes beforehand. Since the physics should be unaffected by the choice of co-ordinates 
to identify points in space, we need to remove this artificially induced anisotropy by averaging over a complete cycle. 
Second, it is unusual to note that the spin- flip parts of (59) are dependent on the x- and y-components of (34), in 
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contrast to the z-component for the Dirac neutrino. The fact that the two remaining components decouple so cleanly 
from the integration indicates a reflection of some fundamental property of the Majorana wave packet when it interacts 
gravitationally. Most interestingly is the fact that a spin-flip term still contributes to the Majorana matrix element 
when ft — > 0, which is not true for the Dirac counterpart. This again follows naturally from the self-conjugation 
property of Majorana neutrinos, though it appears counterintuitive for such terms to appear as written. 



V. GRAVITATIONAL EFFECTS ON NEUTRINO OSCILLATION LENGTH AND ABSOLUTE 

NEUTRINO MASSES 

In order to determine the observable consequences of (56) and (59), we need to obtain the perturbation energy 
which follows from these matrix elements. To demonstrate this explicitly, we make use of the Brillouin-Wigner (BW) 
method [33] of time- independent perturbation theory commonly used in the context of condensed matter physics, as 
opposed to the more well-known Rayleigh-Schrodinger (RS) method [34] we adopted in the previous investigation of 
this problem [12]. In general form, the perturbed energy eigenvalue following the BW method is 

P - -4- M TT l«\ 4- V" IH g mt. |n)| 2 

o + ("I Hint. \n) + 2^ — Z(^)— 

+ ^ ("I #int. \mi) (mi| H iDt . \m 2 ) (m 2 \H iDt . \n) ^ ^ 

m 1 ,m 2 ^n (e - E^™ 1 ^ (^E - E^™ 2 ^ 

where H mt . is the interaction Hamiltonian, \n) is the unperturbed eigenstate, and E^ is its associated energy 
eigenvalue. Although the perturbed eigenvalues from the BW method are generally less accessible than from the 
RS method, given the presence of E in the denominator for the second- and higher-order terms in the perturbation 
expansion, for the special case of a second-order expansion for a two-level system, the perturbed energy eigenvalues 
can be obtained exactly, which is not possible to achieve with the RS method. For this investigation, having a way to 
precisely determine the mass-induced energy difference for a two-neutrino system is extremely useful for accurately 
predicting the effects of spin-gravity coupling on the neutrino oscillation length. 

To determine the corrections to the known neutrino oscillation length due to gravitation, we need to first obtain the 
unperturbed energy eigenvalue E^ from the zcroth-order part of the Hamiltonian Hq. Given (11), it follows that 



where 



H \Mr)) = E^Mr)), (63) 



E^> « y/(hko) 2 + m? (hk Q ) + - —tL z ±- 

r 5 r 6 \ 2 



(hko) 





1 2 




¥ 2 m \ 







and L z is the orbital angular momentum along the gravitational source's axis of rotation. 

It is clear from (64) that the unperturbed part of the Dirac Hamiltonian cannot generate any gravitational cor- 
rections to the neutrino oscillation length, since the mass-dependent part of E^ is not coupled to the gravitational 
potential. However, the situation is different with the introduction of the interaction Hamiltonian _ff$ G given by (39). 
Following (62) and introducing for convenience (with a slight abuse of notation) the unperturbed mass eigenstates 

f (K) + M)± (Dirac), 
l±> = (65) 
I (K) + K))± (Majorana), 

it is shown that the second-order perturbed energy E^ is 



(66) 
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As stated earlier, this expression can be solved exactly for E^\ leading to 

e { £ ] = 1{[(e^ + eP) + (±\h* g \±) 

+ y [(4 ±} + E™) + (±| H* G |±)] 2 - 4 [(4 ±} + (±| H* G |±)) 4 T) - | (T| H* G |±)| 2 ] | , 
which reduces to E^ in the limit as H^ G — > 0. 



(67) 



A. Dirac Neutrinos 



For the case of Dirac neutrinos, the matrix elements to evaluate (67) are 



M 



(±|fl* G |±> - (hk )— [Co + Cim + Czm 2 ] , 
(T\H^ G \±) = (hk ) — sin 9 [D + D im + D 2 m 2 ] 



(68a) 
(68b) 



Upon substitution into (67) and after performing a power series expansion with respect to m, it is shown that the 
perturbed energy is 



4?> = ^™-o) 



2M\ 1 , 

1 + - m 2 

r ) 2 



(ft fco) [(F Dirac ) + (Ff m + (F 2 Dirac ) to 2 ] , (69) 



where 



^Dirac 



Ft 



Dirac 



^iDirac 



^Dirac 



Dirac 
X± 



Co 



M\ 1 /A 



r / 5 V fco r 



C 2 



M\ 1 



5 (fc r) 
5(fc r) 



^ ) - 5 (fco r) 

2 



3/2 



r J \ r 2 

(C 2 + 2C C 2 ) f-y+4sin 2 0(^ 2 + 2 J D o ^ 2 ) " 



\ r J 



(A Di rac) 2 = 25(fc r) 2 C 2 ^ 2 +[l6+100(fc r) 2 sin 2 ^ 2 ] 



(MVLR 2 
\ r 2 



± 40 (fc r) Co 



(70a) 
(70b) 

(70c) 
(70d) 



(70e) 



(70f) 



This leads [8, 35] to the final expression for the neutrino oscillation length L osc . = 2tt / \ E^ — E^J , where 



E (±) _ E (±) 

J - / m2 mi 



(ft fco) 



(^ Dirac ) (ma - mi) + (V 2 Dirac + ±) (to 2 - to 2 ) 



(71) 



It is worth emphasizing again that (71) is formally an exact expression, up to second order in to, for the mass- 
dependent perturbation energy difference, and all plots presented for the Dirac neutrino case are based on (69)-(70f). 
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Nonetheless, given that both M/r <C 1 and MQR 2 /r 2 <C 1 for all numerical analyses considered, it is worthwhile to 
determine the leading order contributions to the mass-dependent energy difference. Therefore, it is shown that 



(±) 



E 



where 



-— — L* + (hk Q )[(l 
+ (»*») ^ t« irac ) 



2M 
r 



1 -2 

-m 



(a**) ^ [(g^t) + ( G ?(±T) ™ + ( G °(±T) ™ 2 



( J Fr i Dirac ) m + (Kf 



(Hk 
+ 



1 



m 9 — m 



M 
r 



( /^Dirac \ 
V G 1(±) J 



+ 



MflR 2 



(m 2 - mi) + 



-™ ^ -iDirac\ 



( r Vi 



(72) 



MflR 2 



(£T 2 Dirac ) 



(m| - mf) 



(73) 



y^rDirac 
°0(±) 



G i(±) 



y^rDirac 
<*2(±) 



y^Dirac 
A 



j^-Dirac 



^Dirac 



(A Dirac ) 2 = 4 + 25(fc 2 sin 2 ^ 2 . 



V ± 2 



5*4 



Co, 



fc r 



Ci - 50 ( ) sin 2 6» Co £>o £>i 



G> - y sin 2 9 [C (D 2 + 2D D 2 )+2 & D £>i 



71(^) sin 2 0Co^ 2 i? 2 



= 5 



5 \k r 
A 



sin 2 0A)-Di, 



5 / fco r 



sin 2 6» 



(£ 2 + .Do D 2 ) 25 2 sin 2 D 2 a D\ 



(74a) 
(74b) 



(74c) 

(74d) 
(74e) 

(741) 
(74g) 



The explicit expression lor the leading-order gravitational corrections to the mass-induced energy difference (73) 
can be grouped into two categories, namely the special case where 6 = 0, it, followed by observations off the axis of 
symmetry where 8 ^ 0, n. For the first case, it is readily shown that 



irac 
(±) 



/-vDirac 
G 2(±) 



= ± 2 C ^ 



(1 ± 1) C 2 



^Dirac j^Dirac 



o. 



= 0, 7T 



(75a) 

(75b) 
(75c) 



19 



and 



x 2 



1 



5 

32 



H |9, 

r 



-2 [CiD -2CbI>i] f 
^ [D? - .Do £> 2 ] + ^ [C 2 D - & Di] 



r 



5 /r/i? 
256 Vr/i?- 



g 3 (67r (J -l)ln(^) 



15 
64 



7r q 



(79a) 



(79b) 



B. Majorana Neutrinos 

When performing the same calculation for Majorana neutrinos, some care has to be taken for evaluating the spin-flip 
part of the matrix element, since it has an explicit dependence on tp. That is, 



<T|tf* G |±> = 

± (hko) sm9 sin^|^(±|cr|T) S [C 0y + C ly rh + C 2y fh 2 ] + M ^ <=bjcr|=F> a [D 0£ + D 1& fh + D 2& m 2 ] j 
= ±(hk ) sin 9 sin ip < cos 9 



sinip ( [C ji + Ciy m + C 2y m 2 ] + cos tp ) [-Dos + D 1$ m + D 2S . m 2 ] 



± i 



cos I — I [C y + C ly m + C 2y m \ - sin tp I — — I [D x + £>ix m + D 2£ m \ 



(80) 



Because the Lense-Thirring metric is axisymmetric, the ^-dependence in (80) needs to be averaged over a complete 
rotation to remove this unphysical anisotropy in the azimuthal direction. Therefore, by setting 



kti^ g i±>i 2 - ^^ 2,r |(TiH* G i±>r 



dp, 



(81) 



it follows that 



(T|#* G |±)| 2 - \ {hk ) 2 sin 2 ' 



^1 - ^ Sin 2 [C0y +Clyfh + C 2 y 

( 1 2 „\ /MOi? 2 \ 2 r 21 
+ ( 1 - — sin 2 6M ( — — \ [Do & + D li; fh + D 2S; m 2 ] 



rn 



(82) 



Following the same procedure as performed for the Dirac neutrino, the perturbed energy eigenvalue for the Majorana 
neutrino is 



_ (±) 4 MOR 2 - 

e£> = ■= — 3~L Z + (hk ) 

5 r 6 



2M\ 1 
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This leads to the corresponding oscillation length expression for Majorana neutrinos, where 
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Similarly, the leading order representation of (83) is 
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As shown for the Dirac case, we can obtain the leading order contributions to (87). Again, we treat as a special 
case the calculations where 9 = 0, n, which leads to 
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C. Analytic Comparison of Dirac and Majorana Neutrino Interactions due to Spin-Gravity Coupling 

A complete numerical analysis of the main results presented here is contained in the next section. However, even 
a preliminary formal analysis shows clear evidence that Dirac and Majorana neutrinos are distinguishable due to the 
presence of spin-gravity coupling. Focussing on (75) and (89) for now, we notice that when the neutrino beam is 
directed along the axis of symmetry it follows that K^" &c and kJ 1 ^' are identically zero, implying that the helicity 
states are unaffected along 6 = 0,7r. For neutrinos with negative initial helicity, both G!? lrac and G!f aj ' are also 
identically zero, while there remains a non-zero contribution of G^ u&c = —Ci/2 which does not exist for G^ a ^'. 

When considering neutrino beams off the axis of symmetry, it becomes obvious from (77) and (90) that there 
are significant differences in ^-dependence which appear, along with relatively trivial differences in the coefficients. 
Furthermore, it follows that K^ 1 ™ are proportional to Dj, which in turn have an overall dependence of (r/R) 2 from 
(78a) and (78b). In contrast, jf^aj. are proportional to Dj$, which arc shown by (91a) and (91b) to be dependent 
on (r/Rf. 



VI. NUMERICAL ANALYSIS 



Having obtained the final analytic expressions for (71) and (85) for the Dirac and Majorana neutrinos, respectively, 
we proceed to express them numerically for astrophysically relevant situations. This provides an opportunity to 
determine the likelihood of observing gravitational corrections to the neutrino oscillation length. One obvious example 
is to perform this analysis for neutrinos emitted from the Sun, while another is to use SN1987A as a test case for 
comparison involving a known supernova source. These examples are considered separately below. Although it is 
not obvious that these examples suggest any possibility of observing these effects, they offer interesting theoretical 
insights about the physical behaviour of Dirac and Majorana neutrinos in a gravitational field. 



23 



A. General Properties 

Before going into detail about the examples considered below, it is worthwhile to discuss some of the general 
properties of the plots common in both cases. All the plots presented below are generated based on the expressions 
found in Appendix C. While the specific examples are, in principle, dependent on the value of the mean momentum, 
it turns out that the plots are insensitive to any large variations of the choice of fco considered. With this realized, 
we present the plots assuming hk ~ 1 MeV throughout. In addition, while there is a formal difference between the 
functions FJ and for each of the cases considered due to the initial negative and positive helicity of the neutrino, 
respectively, the actual numerical differences do not register on any of these plots, where at best we have a one part 
in 10 40 difference due to their initial helicity states. 

B. Gravitational Effects on Neutrinos from the Sun 

We begin with an analysis involving the Sun as the gravitational source. Given that [14] in geometric units 
M = M© = 1.48 x 10 5 cm, R = 6.95 x 10 10 cm, and r = 214.8 R = 1 A.U., it follows that M/r = 9.95 x 1(T 14 for the 
Sun. With the sidereal period of the Sun [36] equal to about 25.38 days, it also follows that MflR 2 /r 2 = 3.07 x lCT 16 . 

The expressions for dimcnsionlcss functions Cj and Cj y arc listed in Figure 2. We recall that Cj contribute to the 
spin diagonal part of the matrix element coupled to M/r, while Cj y for the Majorana neutrino contributes a spin-flip 
part that is also coupled to M/r. For C and C ay in Figures 2(a) and 2(b), respectively, they each have a single peak 
located around q = 1, but are of opposite sign and Cq v is two orders of magnitude smaller in amplitude than Cq. 
The Ci and C\ y in Figures 2(c) and 2(d) behave in a similar fashion, with two exceptions in that they differ by only 
one order of magnitude, and that C\ has a smaller peak that is positive-valued around q = 2. As for C 2 and C 2y in 
Figures 2(e) and 2(f), they each have some non-trivial structure to them. Both functions have main peaks which are 
positive-valued at around 5 x 1CT 1 , where C 2 is about two orders of magnitude larger than C 2y . However, Figure 2(e) 
also has a smaller negative peak at q = 2 and an even smaller positive peak at around q = 3, while Figure 2(f) has a 
slight shoulder at q = 3. All functions Cj and Cj y rapidly decay to zero for q < 10~ 2 and q > 5. 

For the dimensionless functions Dj and Djs with respect to q, they are presented in Figure 3. The first observation 
of note is that their magnitudes are about two orders of magnitude larger than their counterparts for Cj and Cj y . It 
is especially interesting to observe that, contrary to the plots in Figure 2, the Dj$ functions are one to two orders 
of magnitude larger than the Dj. As well, Figures 3(d) and 3(f) which describe D 1S . and D 2 x, respectively, have 
additional structure compared to Di and D 2 with the presence of a negative- valued peak at around q = 2. Like the 
plots shown Figure 2, these plots are non-zero only in the range 10~ 2 < q < 5. 

Using Figures 2 and 3, we can obtain numerical expressions for i^ Dlrac and i^ Maj - with the Sun as the gravitational 
source. We begin with a comparison of F\ for Dirac and Majorana neutrinos presented in Figure 4, where Figures 4(a) 
and 4(b) describe F®" ac and F™*' 1 ' , respectively, for ft ^ 0. It is evident from these first two plots that the gravitational 
field can distinguish between Dirac and Majorana neutrinos, given the distinctive properties of their respective profiles. 
From Figure 4(a), the plot is almost completely negative- valued with a single peak for each choice of 9 in the range 
of 10 _1 < q < 1, while Figure 4(b) shows both a positive- valued peak in the same range of q and a negative- valued 
peak around q = 2. The magnitude of the plots for both cases are on the order of 1CT 11 . In addition, we have for 
further comparison Figure 4(c), which describes F x Maj ' for il = 0. This has the interesting property of a non-trivial 
profile with a single peak around q = 1 for different choices of 9 and a magnitude on the order of 10 - . In contrast, 
the magnitude of irpirac f or — only reaches a maximum of 1CP 20 , which suggests that the rotational part of the 
metric makes the dominant contribution to the predicted gravitational effects on the neutrino oscillation length. 

A similar comparison of F 2 for Dirac and Majorana neutrinos, as shown in Figure 5, is also interesting for a different 
set of reasons. From Figure 5(a), it is evident that the magnitude of F 2 Dlrac decreases as 9 increases, which suggests 
that the rotational part of the metric dampens the gravitational corrections to the quadratic mass term in the energy 
difference due to neutrino mass. This effect does not appear as such for F 2 Maj ' in Fi gure 5(b). However, another 
anomalous effect is present in the plots for small q behaviour, since the function does not vanish as q — ► 0, but appears 
to asymptotically approach a finite and non-zero value on the order of 10 -11 . Although it is difficult to precisely 
determine the cause of this behaviour, it does seem that the presence of rotation in the metric has a role in this effect. 
Comparison with Figure 5(c), which describes F^^' for $7 = 0, seems to confirm this interpretation, as the function 
rapidly vanishes for q < 10 _1 with a magnitude on the order of 10 -14 . Again, this is consistent with the results of 
Figure 4(a), where the magnitude of F 2 Dlrac for Q = is also around 10 -20 . 
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q fdimensionless) 
(a) Co (Sun) 




q (dimensionless) 

(b) C 0il (Sun) 
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(c) Ci (Sun) 




q (dimensionless) 
(d) dy (Sun) 




q fdimensionless) 
(e) C 2 (Sun) 




q (dimensionless) 

(f) C 2y (Sun) 



FIG. 2: Plots of dimensionless functions Cj and Cj y with respect to 5, where the Sun is the gravitational source. While both 
Figures 2(a) and 2(b) peak at around q = 1, they are of opposite sign and Co y is roughly two orders of magnitude smaller than 
Co. Figures 2(c) and 2(d) also have opposite sign and peak at around 5=1, though their widths extend from 10~ 2 to 1, and 
Ci y is only one order of magnitude smaller than Ci. Both Figures 2(e) and 2(f) have positive- valued peaks around 5=1 whose 
widths range from 10~ 2 to 1. However, C2 displays a smaller peak of opposite sign not present in C-z y . 
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(a) D (Sun) 



aj 3.0x10 " 




q (dimensionless) 
(b) A,* (Sun) 




q fdimensionless) q (dimensionless) 



(c) Di (Sun) (d) D m (Sun) 




FIG. 3: Plots of dimensionless functions Dj and Dj± with respect to q where the Sun is the gravitational source. In contrast to 
Co and Co y of Figure 2, Do* is an order of magnitude larger than Do, as shown in Figures 3(a) and 3(b). Similarly, Figures 3(c) 
and 3(d) show that Di* dominates over D\, this time by two orders of magnitude. This is also true for Figures 3(e) and 3(f). 



2G 



5.0x10 



0.0 



























— e = o.Stt 














-- S = OAw 














e = o.3jt 














- — 0.2tt 














— 6 -0.1k 






-— 








— = 





























































Uf 1 10" 2 10"' 10° m 1 



q fdimensionless) 
(a) Ff" Iac (Sun) 



-2 i.oxio - 




g (dimensionless) 

(b) Ff aj ' (Sun) 




FIG. 4: Comparison of fi as a function of q for Dirac and Majorana neutrinos assuming various angles of the neutrino beam 
off the axis of symmetry, where 9 is defined by Figure 1. For both Figures 4(a) and 4(b), the plot profiles show that magnitude 
of F\ is most pronounced for 6 = tt/2. For a non-rotating star (fi = 0), Figure 4(c) shows that there still exists a non-trivial 
contribution to i^ 1 Ma j' i though reduced by three orders of magnitude compared to Figure 4(b). 
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FIG. 5: Comparison of F2 as a function of q for Dirac and Majorana neutrinos for varying neutrino beam orientation. Unlike 
for Figure 4, an increase of 6 on Figure 5(a) leads to a decrease in the magnitude of F^** 2 , suggesting that rotation of the 
gravitational source serves to dampen the quadratic mass correction in the energy difference. A more intriguing effect appears 
in the small q behaviour found in Figure 5(b), which does not dampen to zero as q — > 0. This is somewhat surprising, since 
the profile of the matrix element components which contribute to this plot for F 2 aj ' all vanish for small q. For the case of a 
non-rotating star (SI = 0), Figure 5(c) yields a non-zero plot, again with a reduction of three orders of magnitude compared to 
Figure 5(b). 
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Predicted Effects on Neutrinos from SN1987A 



Using SN1987A as the gravitational source [37], M 
: 1.5 1 2 x 10 23 cm, leading to M/r = 1.3704 x 10 ' 



-18 



1.4M = 2.072 x 10 5 cm, R = 10 km, and r = 49 kpc 
Unlike the Sun, SN1987A has a much higher rotational 



-37 



frequency, with an estimated period of 2.14 ms. This results in MflR 2 /r 2 = 8.8702 x 10 

Plots of the functions Cj and Cj$ are found in Figure 6. Essentially, they possess the same properties as found in 
Figure 2, except for an increase in their magnitudes. The corresponding plots of Dj arc found in Figure 7, and they 
also behave in the same fashion as found in Figure 3. As for Dj£, they need to be treated with greater care. This 
is because, as stated earlier, Dj$ are proportional to (r/R) 3 as opposed to (r/R) 2 for Dj (see Appendix C), which 
implies that the expression (MQR 2 /r 2 )Dj£ is proportional to r/R. For r>fl, the contribution oiDj± overwhelms all 
other contributions to the matrix element. This appears to be physically unrealistic, since the effect of rotation should 
decrease with increasing radial distance from the source. Therefore, to compensate for this anomaly, we introduce by 
hand an overall regulator (C.l) to account for the effects of wave packet spread, such that 
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(93) 



which amounts to exponentially damping these functions to zero for large r. 

The plots of i?\ Dirac and F™**' for varying 9 are listed in Figure 8. Like that shown with Figure 4 involving the Sun, 
Figure 8(a) shows that the magnitude of F 1 Dlrac grows with increasing 9 and is on the order of 10~ 2 , while Figure 8(b) 
has a magnitude of order 10~ 4 . Similarly, the comparison of F 1 Dlrac and F^ 1 ^' as a function of shown in Figure 9 
resembles that of Figures 5(a) and 5(c), with magnitudes of 0.5 — 0.6 and 8.0 x 10~ 4 , respectively. 



D. Observational Possibilities 



Because (71) and (85) have both a dependence on Am^i 



m 2 and A77121 = TO2 — mi, where by convention 



we assume m,2 > mi, it follows that we have enough information to algebraically solve for the absolute masses mi 
and m-2 by a parameter fit of q, Am 2 i, and Am|i to any data which registers gravitational corrections to the neutrino 
oscillation length. This leads to 
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It is interesting to note that, while any precise determination of the linear mass difference must come from observation, 
it is possible to determine a theoretical upper bound to Am,2i strictly from consistency arguments that follow from 
assuming that mj > 0. Since it must be true from (94) that 
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(96) 



it follows that we can define an experimental parameter e as some deviation from the presumed mass-induced energy 
difference due to special relativity, such that 
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where AE = E^ — E^' . Then upon solving for e, leading to 
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FIG. 6: Plots of dimensionless functions Cj and Cj y with respect to q, where SN1987A is the gravitational source. Apart from 
their magnitudes, Figures 6(a)-6(f) behave in the same general way as their counterparts in Figure 2. 
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FIG. 7: Dimensionless functions Dj as a function of q for the case of SN1987A. When the effects of wave packet spreading 
are taken into account, it happens that the corresponding Dj^ are exponentially damped to zero, and therefore make no 
contribution to the predicted neutrino energy difference. 
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FIG. 8: Comparison of Fi as a function of q due to the SN1987A gravitational source, for varying neutrino beam angle i 
the behaviour shown in Figure 4. the effect of increasing 6 is to increase the mag nitude of F\ Dirac and i*\ Maj '. 
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FIG. 9: Comparison of F2 as a function of q for Dirac and Majorana neutrinos with SN1987A as the gravitational source. 
While less pronounced as that demonstrated in Figure 5, the effect of increasing 8 on Figure 5(a) also leads to a decrease in the 
magnitude of j , ° lrac ) with a shift of the peak to lower values of q. Unlike that of Figure 5(b), however, the small q behaviour 
of Figure 9(b) is as expected. 



it follows that use of (96) leads to the upper bound for the absolute neutrino mass difference, 

Fi 



Am 2 i < 



(99) 



Figures 10 and 11 describe (99) for the solar and SN1987A cases, respectively, as a function of q and for a variety 
of choices for e. We notice immediately that the upper bound for A?fi2i is essentially constant in q for most choices 
of e, except when q > 10. Direct comparison of Am!? 1 lrac with Afh^^' shows that the former allows for a much 
higher upper bound with even very tiny choices of e, while the latter is very sensitive to the choice of e and leads 
to more restrictive upper bounds for Am2i- For the case of the Sun, Figure 10(a) shows that Am 2 3 1 lrac < 0.491 for 
10 -30 < e < 10 -16 , which then increases to Am^ 1 ™ < 0.54 for most choices of q when e = 10 -14 . Given that 
the mean momentum Hko for solar neutrinos is on the MeV scale and we assume oscillations between electron and 
muon neutrinos, the upper bound for Am^ mc = Am^'^fWo) is largely consistent with the upper bound of m v < 
0.170 MeV [7], as determined by experiment. In contrast, Figure 10(b) shows that Afn^^' < 10~ 3 for e = 10~ 10 
and Am 21 a ^' < 10~ n for e = 10 -2 , while Figure 10(c) leads to Am" 3 '' < 10~ 6 and Am 21 aj ' < 10 -14 for the same 
respective choices of e. These results suggest the prediction that Majorana neutrino masses must be closely spaced 
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FIG. 10: Upper bound of Arfi2i for various choices of e, where the Sun is the gravitational source. It appears that Arn,2i is 
largely a constant in q throughout the whole range available. Comparison of Figure 10(a) with Figure 10(b) shows that the 
mass separation for Dirac neutrinos can be several orders of magnitude larger compared to that of Majorana neutrinos. From 
Figure 10(c), the absence of £1 serves to reduce the upper bound of Aifi2i for each choice e by three orders of magnitude. 
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FIG. 11: Upper bound of Arn.21 as a function of q for the case of SN1987A. From Figure 11(a), it appears that the maximum 
mass separation for Dirac neutrinos is independent of the choice of e, while Figure 11(b) shows that the upper bound decreases 
in proportion to each order of magnitude increase of e. 



together to be self consistent with this model, while Dirac neutrinos have the option of also having a much larger 
mass separation. This property is also reflected in Figures 11(a) and 11(b), where Afh^" ac < 0.074 for the full range 
of q considered, while Ara^ a '' < 10~ 6 for e = 10~ 3 and reduces to Ami^' < 10~ 8 for e = 10 _1 . Furthermore, a 
much larger e is required to obtain Am 21 a ^' than for Aml? 1 lrac , which immediately follows from the fact that _F 1 Dlrac 
carries opposite sign compared to i 7 f )lrac , while F^ 1 ^' and F^ 1 ^' have the same sign. 

Finally, we consider the likelihood of observing the gravitational corrections in future neutrino oscillation experi- 
ments. The best possibility appears to be with the efforts of SNO [38] and Borexino [39] to measure low-energy solar 
neutrinos on the order of 1 MeV, whose flux rate is many orders of magnitude larger than for previous measurements. 
This can be done in combination with a precision measurement of ArnJi from an independent observation, such as 
from the T2K [40] long baseline neutrino oscillation experiment under development. In particular, SNO is focussing 
attention on measuring the monoenergetic solar neutrinos [41] with energy E v — 1.44 MeV due to the pep reaction 
component in the pp — I chain [7], with an uncertainty of 3-5% [41]. The advantage with this measurement comes 
from knowing that this reaction is confined to a very narrow energy range and suggests a realistic value for q that 
falls within the non-zero values for Fj. As for neutrinos from supernovae like SN1987A, while the plots from Figures 8 
and 9 show a significant spin-gravity correction to the mass-dependent energy difference, it seems unlikely [41] that 
we may be able to find suitable candidates which can provide the required precision to identify a recognizable signal. 
However, it remains to be seen whether that is the final answer. 



VII. CONCLUSION 



We have shown in this paper that spin-gravity interaction, via gravitational phase and time-independent pertur- 
bation theory, offers the possibility of distinguishing between Dirac and Majorana neutrinos propagating as wave 
packets. The distinctions between the two neutrino types are unambiguous at the matrix element level, and lead to 
meaningful predictions in gravitational corrections to the neutrino oscillation length. This has been explicitly demon- 
strated using both the Sun and SN1987A as test cases for detailed analysis. The results presented here demonstrate 
that gravity can be an interesting probe to better understand the fundamental nature of neutrinos and possibly other 
subatomic particles. Because the associated matrix elements for the two neutrino types generates a linear mass term 
in the perturbed total energy that is coupled to the background gravitational field, we have the theoretical possibility 
of determining the absolute masses for a two-flavour oscillating system. Furthermore, this model can be used to 
determine for both Dirac and Majorana neutrinos the allowable upper bound for the absolute mass difference, and 
shows that while Dirac neutrinos can have a large mass separation, Majorana neutrinos must necessarily have a much 
smaller mass separation under the same conditions. 

It is worthwhile to identify details where improvements of the model can be employed. For instance, we made 
the assumption that the width of the neutrino wave packet is constant in time. Since this is strictly not true, it 
will be necessary to incorporate the spreading of the wave packet in any future development of this project. In so 
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doing, however, we change the fundamental nature of the problem such that we then require either time-dependent 
perturbation theory or some other compatible technique. Another technical point to address concerns the use of a 
regulator to tame some isolated terms in the matrix elements which otherwise grow infinitely large in the spatial 
integration. This issue may actually be resolved if a dynamical wave packet width is accounted for in a more refined 
model. Since we know of at least three neutrino flavours in existence, further development of this project may include 
a three-flavour oscillation, which introduces considerable complexity to the original problem. On a more fundamental 
level, it is possible to consider abandoning the wave packet model altogether in favour of a quantum field theory 
approach in curved space-time [42] . It remains to be seen what is the best route to follow for any future developments 
of this analysis. Nonetheless, it seems clear that the results presented here merit further investigation. 
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APPENDIX A: EVALUATION OF MOMENTUM INTEGRALS 

The required integrations over k are described by the functions 
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where 
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and n are integers and —2 < n < 2. For q = ko/a p , the evaluation of (A.3) for non-negative n is relatively 
straightforward, the results of which are 
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In order to allow for the final integration over r, the complex error functions present in (A.4)-(A.6) and elsewhere 
have to be approximated for suitable limiting cases. For a complex variable z = x + iy with \z\ S> 1, it is generally 
true that 
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If x = q/2 and y = k$r / q, then for all the choices of q considered for given kg and R, it is always true that y 1, 
and that x/y = | (q 2 /(k^r)) <§C 1. By performing an expansion of (A. 7) with respect to x/y, it can be shown that to 
leading order for fixed k , 
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Substitution of (A.8) into (A.4)-(A.6) leads to (A.l) and (A. 2) once the real and imaginary parts of the final expressions 
are evaluated. 

The integrals Q{— 1) and Q(— 2) are more challenging to perform, given that k is in the denominator, suggesting 
the possibility of infrared divergences for these expressions. However, there is a procedure available to obtain these 
integrals without any troubling behaviour as k —* 0. To accomplish this, we first define a function f„(ko) where 
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such that Q c {—n) = Re [/„(fc )] and Q s (—n) = Im [/„(fc )]. By differentiating (A.9) with respect to fco, it follows that 
the integral must satisfy the first-order differential equation 
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which immediately leads to the solution 
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with constant of integration K n . Since /„(fc ) — > as fc — > oo, it is clear that K n — from this boundary condition. 

By repeated use of (A. 11), f n {k ) can be explicitly evaluated for all n > 1, provided that an explicit solution for 
/l(fco) is obtainable. Fortunately, given that /o(fco) = <3(0) by construction, /i(fco) can be expressed upon substitution 
of (A.4) into (A. 11). Utilizing (A.8) and the integral 
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Finding (A. 11) for n = 2 is also a straightforward process, given (A. 14). This leads to 
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Therefore, the final expressions for Q c (n) and Q s (n) are the following: 
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APPENDIX B: EVALUATION OF RADIAL INTEGRALS 



For this paper, most of the radial integrals have the form 
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where a is a positive integer and is an integer in the range —a < < a. For nonnegative n, (B.l) can be immediately 
integrated, leading to 



S(a,0,O) 
S(a,0,l) 

S(a,0,2) 



2 \ a Uo 



e "/JV/(4«) erf 



kor\ . fq 



kn I 2a 



JL J _ e »/3feor e -a(fc 0j R/<?)" - {til. / _ e -r?7(4«) er f 



9 



/~7T 



(1) 



2a \ 5 / 4a 2 



4 V a 



e i0k o r e -a(k R/q) 2 



1 JL(2a-0 2 a 2 )e-^ 2 /^ erf 
8 V a 



(to 

V 9 



r\ . fq\ 



^a 



(B.2) 



(B.3) 



(B.4) 



where 

It then follows that 



erf 



fq 



(I) 



1 f 1 



la 7r V fco r 



1 + i 



0q ( q 



2a V fco r 



e -/3V/(4a) e -a(k R/qf £ ipkar _ 



S(a,/?,0) 



2a V^o-R 



cos(/3 fcoi?) 



0q ( q 



2a V fcoi? 



sin(/3 fc i?) 



+ i 



£ ( fcS?) cos(/3 fi) + sin(/? fcoi?) } e_Q(feoK/9)2 ^ 



(B.6) 



S(a,0,l) = 



2a V^oi? / 



5v xy 

4a 2 I fc i?/ 



+ i 



g^cos^fcoi?) 



2 q 2 ( q V 
4a 2 ^fc i?y 



sin(/3fc i?) 



sin(/3 fc i?) 

e -a(feoi?/<?) 2 R 2 



(B.7) 



S(a,/?,2) 



2a \k Rj 



-0q 




2a 


\k RJ 



cos(0 k R) + sin(0 fcoi?) 



e -a(k R/q) 2 R 3 



(B.8 



However, for positive-valued n, the evaluation of (B.l) is less straightforward. Fortunately, by treating as a 
continuous parameter and making use of the identity 
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it is possible to solve for S(a, 0, —1) using (B.6). This leads to the iteration solution 
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As expected, when /3 — > — /3 the real parts of (B.6)-(B.8) and (B.10) remain unaffected, while the imaginary parts 
change sign under this transformation. 
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APPENDIX C: EVALUATION OF THE DIMENSIONLESS FUNCTIONS Cj AND D 3 



In this Appendix, we list in analytical form the dimensionless functions Cj and Dj for the Dirac and Majorana 
matrix elements, due to the gravitational phase <£>g- Up to the approximation of the complex error functions, as 
described in Appendices A and B, these expressions for Cj and Dj are exact. While almost all the terms in the final 
expressions below are well behaved for the full range of q, there are some isolated terms which diverge when q — > 
and/or when r — > oo. For these terms, we insert by hand an exponential decay regulator of the form 



£(q,k ,r,R) = exp 



k%{r-Rf 



(C.l) 



which mimics the effects of wave packet spreading with time to provide a physically plausible cutoff of the respective 
integrals. The choice of (C.l) results from evaluating the time evolution of a one-dimensional wave packet with 
ultrarclativistic energy, following standard techniques [11]. 
The Cj and Dj can be written in the general form 
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Dj (q, k , r, R) = {Dj) ab cos a (ko R) sin 6 (ko R) , 



(C.2b) 



where (Cj) b and (Dj) b are expressed in terms of q, the dimensionless parameter /j, = qj (ko R) <C 1, and (C.l). We 
proceed to list the components of (C.2) for the Dirac and Majorana matrix elements, respectively. 



1. Dirac Neutrinos 



For the Cj functions, we have the respective components (Cj) ab , where 
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For the Dj functions, the components (Dj) ab are 
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2. Majorana Neutrinos 

The dimensionless functions Cj y and Dj$ which define the Majorana matrix element can be expressed as shown in 
(C.2). Therefore, it follows that the components for Cj y are 
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